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An extensive debate on quantum non-demolition (QND) measurement, reviewed in Grangier et 
al. [Nature, 396, 537 (1998)], finds that true QND measurements must have both non-classical 

a* 

state-preparation capability and non-classical information-damage tradeoff. Existing figures of 
merit for these non-classicality criteria require direct measurement of the signal variable and are 
qq ■ thus difficult to apply to optically-probed material systems. Here we describe a method to demon 

in . 

\Q , strate both criteria without need for to direct signal measurements. Using a covariance matrix 

— -^ [ formalism and a general noise model, we compute meter observables for QND measurement triples, 

which suffice to compute all QND figures of merit. The result will allow certified QND measurement 

of atomic spin ensembles using existing techniques. 
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I. INTRODUCTION 



A quantum non-demolition (QND) measurement is one which provides information about 
a quantum variable while leaving it unchanged and accessible for future measurements. 
The approach was originally suggested as a means to avoid measurement back-action in 
gravitational wave detection |ll45|. QND measurements of optical fields both provided the 
first demonstrations and led to a considerable refinement of the understanding of QND 
measurements in practice [6[. More recently, QND measurements have been emp loyed to 
prepare spin-squeezed atomic states [7H11J and with nano- mechanical systems 12 1. 

In a generic QND measurement, a 'meter' and a 'system' variable interact via a selected 
Hamiltonian. The meter can then be directly measured to gain indirect information about 
the system. In the context of optical QND measurements, the question of when a mea- 
surement should be considered QND has been much discussed (see |6( for references). Two 
distinct non-classicality criteria emerge: A state preparation criterion requires small uncer- 
tainty in the system variable after the measurement while a second criterion describes the 
information-damage tradeoff in the measurement. While some operations such as filtering or 
optimal cloning can be non-classical in one or the other criterion, a true QND measurement 
is non-classical in both |6|. 

With the aid of figures of merit 13Ml5| describing the quantum-classical boundary, opti- 
cal QND measurements satisfying both criteria have been demonstrated [l4|, I16M24J . These 



figures of merit make use of the fact that the optical signal beam, after the QND measure- 
ment, can be verified by a direct, i.e., destructive, measurement with quantum-noise-limited 
sensitivity. Typically such a direct measurement is not available in atomic QND. Rather, re- 
peated QND measurement has been used to show the state preparation criterion |8l-ll ll. u 
by conditional variance measurements. Here we show how repeated QND measurements 
can also be used to test the information-damage tradeoff, and thus to certify full QND 
performance without direct access to the system variable. 



II. MODEL 

As in the pioneering work by Kuzmich, et al. |7|, |26[, we consider the collective spin 
of an atomic ensemble, described by the vector angular momentum operator J. We note 



that a variety of other physical situations are described in the same way, e.g. by using a 
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pseudo-spin to describe a clock transition J9| . The optical polarization of any probe pulse is 
described by a vector Stokes operator S 

Si = -aVja, (1) 

i = x,y, z where <Tj are the Pauli matrices, a = {a + , a_} T and a± are annihilation operators 
for circular-plus and circular-minus polarizations. 

We define Stokes operators P, Q for the first and second pulses, respectively. The oper- 
ators J, P, Q each obey the angular momentum commutation relation [L x , L y ] = iL z and 
cyclic permutations (for simplicity, we take h = 1). For notational convenience, we define 
the combined optical variables C = P © Q and the total variable T = J © C. We will be 
interested in the average values of these operators, which we write as J = (J) and similar, 
and the covariance matrices, which we write as 

J = -(3A J + (JAJ) T )-(J) A (J) (2) 



27|,|31 



and similar. Our approach follows that of Madsen and M0lmer 

We assume that the input probe pulses are polarized as _PJ m ) = Q^ = S^ and that the 
other average components are zero. We take the initial covariance matrix for the system to 
be 

T = J © C (3) 

This form of the covariance matrix allows for arbitrary prior correlations (including corre- 
lated technical noise) among the two optical pulses, but no prior correlations between the 
atoms and either optical pulse. 

The interaction is described by an effective Hamiltonian 

H cS = gJ z S z , (4) 



where g is a constant 28j. This QND interaction, to lowest order in gr, where r is the 
interaction time of the pulse and atoms, produces a rotation of the state, T^ out ^ = T^ — 
ir[T^ n \ H e s\. This has the effect of imprinting information about J z on the light without 
changing J z itself: 

S( out ) = S^ n) + k' S<£ n) Jf^ (5) 

jjout) = jjin)_ (Q) 



Here k' = gr and S is P or Q depending on which pulse-atom interaction is being described. 
The rotation can be described by a linear transformation T^ out ) = MpT^ in ^ (and thus T( out ) = 
MpT^Mp) where Mp is equal to the identity matrix, apart from the elements (Mp) 2: Q = 
k! ' J^ n \ an d (Mp) 5i3 = n'S^ n \ For later convenience, we define k = k'S^ = grS^ n \ 
The effect of the second pulse is described by the matrix Mq = XMpX where 

/ I N 
X= 1 ® J 3 (7) 

exchanges the roles of P and Q, and 1% is the 3x3 identity matrix. 

III. REDUCTION OF UNCERTAINTY BY QND MEASUREMENT 

We first consider the case in which the interaction does not introduce additional noise 
(although both the input atomic and optical states may be noisy). After interaction with 
the first pulse, but before the arrival of the second pulse, the state is described by T P = 
MpToMp. A component P y of the first pulse is measured. Formally, this corresponds to 
projection along the axis mp = {0, 0, 0, 0, 1, 0, 0, 0, 0} T , and Tp is reduced to 

f PD = f P - T P (Tl Q fpTl Q ) MP fT = f P - f P U Q f^/Ti[U Q fp} (8) 

where Uq = mp A mp is the projector describing the measurement and () MP indicates the 
Moore-Penrose pseudo-inverse. 

We can directly calculate the resulting variance of J z , 

E[v a r(J z )\P y ] = (T PD ) 3 ,3 = Js,s 2 j° 2 ^ ■ (9) 

« ^3,3 + ^2,2 

This has a natural interpretation: The variance of the detected projection P y has two con- 
tributions: k 2 Js^ from the atomic signal and 6*2,2 from the pre-existing optical noise. 73,3 
is reduced by the factor 1/(1 + SNR) where SNR is the signal-to-noise ratio of the measure- 
ment. A similar result is found in reference J27J. This post-measurement variance of the 
signal variable describes the state-preparation capability of the QND measurement. Absent 
the ability to directly measure J z , we must look for observables which contain this same 
information. 



IV. OBSERVABLE CORRELATIONS 

After interaction with both the first and second pulses, we have Tpq = MqTpMq. This 
matrix contains the variances and correlations that are directly measurable, namely those 
of the two light pulses. These are 

var(P„) = C 2 , 2 + k 2 J 3 , 3 (10) 

var(g s/ ) = £5,5 + /€ 2 J 3 , 3 (11) 

cov(P y ,Q y )=C 2 , 5 + K 2 J 3 , 3 . (12) 

We note that for k — 0, e.g. if the atoms are removed, the values are 

var NA (P y ) = 6*2,2 (13) 

vai NA (Q y ) = C 5i5 (14) 

COV NA (Py,Qy) = C 2 ,5- (15) 

We see that the state preparation capability can be expressed in terms of measurable 
quantities as 

E[^(J z )\P y ] = J 3 /-^^, (16) 

which uses the variance of the two measurements to determine the SNR. Another formula- 
tion, 

E[^(J z )\P y ] = J 3 ,3 ^NA(P y ) ^^ (1?) 

ySLT^A{Py> + COY{Py, Q y ) - COV NA (P y , Q y ) 

expresses the residual variance in terms of the atomic contribution to the correlation between 
first and second pulses. 

These simple expressions are only valid for noise-free interactions, however. In a real 
experiment, other effects are present which introduce both noise and losses in the atomic 
and optical variables. We now account for these other effects. 

V. GENERAL NOISE AND LOSS 

We now consider noise produced in the atom-light interaction itself, as well as losses. The 
noise model we employ is very general. The interaction of the first pulse with the atoms is 
described by 

f P = Mpf Ml + N P (18) 



We assume that the coherent part of the interaction is Mp = r^/3 © tlI 3 © ^3 apart from 
the elements [Mp)i,% = ^xi an d (-Wp)s,3 = ~ K ^ X ■ Here Tai^l describe the fraction of 
atoms and photons, respectively, that remain after the interaction. Thus Mp includes both 
the effect of H e Q and linear losses. We leave Np completely general, except that it does not 
affect Q: Np = N © O/3, where N is a six-by-six symmetric matrix. 
Similarly, we describe interaction with the second pulse as 

f PQ = M Q f P M T Q + N Q (19) 

where M Q = XM P X and N Q = XN P X. 

Note that we assume that both the interaction M and the noise N are the same for the 
first and second pulses (but act on different variables, naturally). This implies that optical 
characteristics of the pulses such as detuning from resonance are the same, a condition that 
can be achieved in experiments. It also assumes that the noise generated by the interaction 
is incoherent and state-independent, as opposed to a more general, state-dependent noise 
iV(J, S). Nevertheless, in many situations J and S are nearly constant (only small quan- 
tum components change appreciably), so that any reasonable N(3, S) would be effectively 
constant. 

As above, we can directly calculate Tpp, and Tpq to find 

E[y & r(J z )\P y ] = J 3;3 r\ + N 3 , 3 - (^.s + iVs,) 2 (2Q) 

k 2 J 3l 3 + riC 2 ,2 + N 5j5 
and 

var(P y ) = r 2 L C 2a + n 2 J 3 , 3 + iV 5 ,5 (21) 

var(Q„) = r 2 L C 5>5 + K 2 (ry 3 ,3 + N 3 , 3 ) + iV 5 ,5 (22) 

cov(P y , Q y ) = r 2 L C 2 ^ + K 2 (r A J 3t3 + N 3>s /k). (23) 
Equation (TT3|) still holds for the case with no atoms. We define 

5var(P y ) = vax(P y ) - vai NA (P y )r 2 L (24) 

5vai(Q y ) = vai(Q y ) - vai NA (Q y )r 2 L (25) 

5COV(Py, Qy) = COv(Py, Q y) ~ COV N A (Py , Qy)^, (26) 

where the r L factors are included to account for atom-induced optical losses. It is then 
simple to check that 

E[var(J M )\P y ] = J 3i3 + k- 2 (sv & r(Q y ) - 5var(P,) - Scav2 (Qv> p v) \ . (27) 

V var lP/J / 



We note that the QND measurement reduces the variance of J z if the quantity in paren- 
theses is negative, i.e., if 

5cov 2 (Q y , P y ) > var(P y )[5vai(Q y ) - <Jvar(P„)]. (28) 

Again, there is an intuitive explanation: 5cov(Q y ,P y ), which arises from the fact that 
both pulses measure the same atomic variable J z , is a measure of the atom-light coupling. 
[5vai(Qy) — 5vai(P y )] expresses the difference in atom-induced noise between the first and 
second pulses. This difference indicates a change in the atomic state, namely an increase in 
var( J z ). The condition of Equation (128]) compares these two effects and can be tested know- 
ing the statistics of the various measurements on S y and the optical transmission tl- The 
factors k 2 , J3 : 3 in equation (}2"T|) must be determined by independent means. For example, k 
can be found by measuring the rotation of a state with known ( J,) ^ and J33 from the 



number of atoms, or the observed noise scaling of a known state 29|, [30] 



VI. THREE-PULSE EXPERIMENTS 



The above description of two-pulse experiments can be extended straightforwardly to 



three or more pulses [31]. While a two-pulse experiment, plus prior knowledge of k and 
J$ t 3, gives sufficient information to find the post-measurement variance, and thus test the 
state-preparation property, a three-pulse experiment is required to find the other quantities 
used to characterize QND measurements. 

If R denotes the Stokes vector of the third probe pulse, then statistics such as var(P y ) and 
cov(P y , Ry) can be determined, and these in turn provide enough constraints to determine 
the loss and noise. Expanding our system toT = J©P©Q©R, and defining interaction 
and noise operators M R , N R in the obvious way, a direct calculation finds several useful 
relations 

= Scov(P yj R y ) 

SCOV(Py,Qy) y ' 

2 = Swsi{Ry) - 5vai(Q y ) 
Ta 5vai(Qy) - Svai{P y ) { ' 



k 2 N 3:3 = 5var(Q y ) - 5var(P y ) + k 2 J 3 , 3 (l - r 2 A ) (31) 

^iV"3,5 = Scov(P y , Q y ) - n 2 J 3t3 r A (32) 

iV 5 , 5 = 5var(P,) - k 2 J 3i3 . (33) 



VII. MEASURES OF QND PERFORMANCE 

To quantify QND performance, Holland et al. use the degree of correlation between 
various combinations of the input and output system variable X = J z and meter variable 



Y = S y variables [13] • They define three figures of merit, each of which is unity for an ideal 
QND measurement. These describe the measurement quality, the preservation of the initial 
value, and the state preparation capability, respectively: 

COV 2 (X in ,y° ut ) K 2 jj fl « 2 J 3 ,3 , , 

x ™ Ymt - var(X-)var(y-t) j^(f P )^ var(P,) l > 

C 2 _ cov 2 (X in ,X out ) _ r\Jl 3 



n ' Xout " var(X^)var(X-t) ~ J 3>3 (f P ) 3i3 

K 2 J^ 3 5 COY 2 (Py.Ry) 



r 2 

Wj^out yout 



8cov 2 (P y ,Q y )[Svai(Qy) - 5vai(P y ) + n 2 J 3 - 3 ] 

cov 2 (X° ut ,y° ut ) (f P )l 5 

var(X- t )var(F- t ) " (f P ) 3 , 3 (f P ) 5 , 5 
5cov 2 (P y ,Q y ) 



(35) 



yar(P y )[5vai(Q y ) — 5v&r(P y ) + k 2 J 



(36) 
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VIII. NON-CLASSICALITY CRITERIA 



Q 



Roch, et al. |l4[ and Grangier et al. 15( define non-classicality criteria using the con- 



ditional variance AX 2 ,, as in Eq. ( 127)) . and the quantities AX^, the measurement noise 
referred to the input and AX 2 , the excess noise introduced into the system variable. All are 
normalized by the intrinsic quantum noise of the system variable, a quantity which may de- 
pend on the system or the application. For example, in a spin-squeezing context the natural 
noise scale is Jo = I (Jr) 1/2 = ^3,3, the J z variance of the input x-polarized coherent spin 
state, i.e., the projection noise. Here we choose to normalize AX^ by J , and AX 2 , m , AX 2 
by Tj\ J , reflecting the reduction in size of the spin due to losses in the measurement process. 
The relation of information gained to damage caused is non-classical if AX s AX m < 1. We 
find 



AX 



2 _ E[vax(J z )\P y ] 

\m j 

taJo 



5cov(P y , Q 



!)) 



SCOV(Py, Ry) 



l + U 2 



Jo)" 1 (sv*r{Q v ) - 5^{P y ) - 6C °l2^ v) ) 



(37) 



_ C 2 ,2r 2 L + ^5,5 _ var(P J/ ) - k 2 J 3 , 3 



AX 2 = ~^-l . -0,0 = —^vj _■» ^ (3g) 

K Z J K A J 

AX 2 _ (fp) 3 ,3 - ^3,3 = 6cov(Py,Q y )[5vai(Qy) - $v&r(Py)] ^ 

r A J 5cov(P y ,R y )n 2 J 



IX. CONCLUSIONS 

Using the covariance matrix formalism and a general noise model, we have shown that full 
certification of QND measurements is possible without direct access to the system variable 
under study. We find that repeated probing of the same system gives statistical informa- 
tion sufficient to quantify both the state preparation capability and the information-damage 
tradeoff. The results enable certification of true quantum non-demolition measurement of 



y 



applicable to ongoing experiments u sing QND measure 



material systems, and are direct 

ments for quantum information 28J and quantum-enhanced metrology 30 
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X. ADDITIONAL MATERIAL 

The calculations described in this article can be performed in Mathematica using the 
notebook "ThreePulseCMCalculator," available as an ancillary file. 
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